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Abstract Star formation is thought to be triggered by 
gravitational collapse of the dense cores of molecular 
clouds. Angular momentum conservation during the 
collapse results in the progressive increase of the cen- 
trifugal force, which eventually halts the inflow of ma- 
terial and leads to the development of a central mass 
surrounded by a disc. In the presence of an angular 
momentum transport mechanism, mass accretion onto 
the central object proceeds through this disc, and it 
is believed that this is how stars typically gain most 
of their mass. However, the mechanisms responsible 
for this transport of angular momentum arc not well 
understood. Although the gravitational field of a com- 
panion star or even gravitational instabilities (particu- 
larly in massive discs) may play a role, the most general 
mechanisms are turbulence viscosity driven by the mag- 
netorotational instability (MM), and outflows acceler- 
ated ccntrifugally from the surfaces of the disc. Both 
processes are powered by the action of magnetic fields 
and are, in turn, likely to strongly affect the struc- 
ture, dynamics, evolutionary path and planet-forming 
capabilities of their host discs. The weak ionisation of 
protostellar discs, however, may prevent the magnetic 
field from effectively coupling to the gas and shear and 
driving these processes. Here I examine the viability 
and properties of these magnetically-driven processes 
in protostellar discs. The results indicate that, despite 
the weak ionisation, the magnetic field is able to couple 
to the gas and shear for fluid conditions thought to be 
satisfied over a wide range of radii in these discs. 
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1 Introduction 

The birth of stars, a process so commonly occurring 
throughout the whole universe, is as puzzling as it is 
fascinating. The current paradigm for star formation, 
framed originally by Kant and Laplace in the 18 th cen- 
tury, suggests that stars are born via gravitational col- 
lapse of the dense cores of molecular clouds, in turn 
condensed out of the diffuse interstellar medium. Un- 
fortunately, this compelling visualisation and its phys- 
ical underpinning are in apparent contradiction. 

The reason for the above statement is that conser- 
vation of angular momentum during the collapse re- 
sults in the progressive increase of the centrifugal force, 
which eventually halts the infalling gas and leads to 
the development of a central mass (i.e. a 'protostar') 
surrounded by a flattened disc of material (an 'accre- 
tion disc'). Observational evidence for the presence of 
such discs around young stellar objects is compelling. 
It comprises imaging in the near infrared and optical 
wavebands (e.g. see the reviews by [Watson et al. |2007| 



McCaughrean et al. 2000) as well as interferometric 



studies that have resolved the velocity profile and struc- 
ture in the inner regions of discs, up to ~ a few tens 



of AU from the centre (e.g. Dutrey et al. [2007 Wilner 
& Lay 2000). This phenomenon is, in fact, not lim- 



ited to young stars. Many astrophysical systems exhibit 
the characteristic disc-like structure that naturally re- 
sults when inward motion in the plane of rotation is re- 
stricted by angular momentum conservation, while col- 
lapse continues in the perpendicular (polar) direction. 
Such structures are commonly associated, for example, 
with the discs of material feeding the cores of active 
galaxies and black holes. 

The difficulty in progressing past this stage is clear 
when we recall that these protostellar discs are dynam- 
ically stable. A typical disc is differentially rotating, 
with an angular momentum (L) profile that increases 
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with radius (dL/dr > 0). Like our own planet, its mate- 
rial is unable to move closer to the central object unless 
it loses angular momentum. This can only be accom- 
plished by redistributing it to other parts of the disc, 
or by transferring it away from the star-disc system al- 
together. Indeed, the efficiency and evolution of this 
phase of the star formation process is regulated by the 
rate at which angular momentum can be transferred 



away from the accreting matter (e.g. Adams & Lin 
1993). Fortunately for us, the Earth is unable to do 



so (and stays in its orbit!) but, how does nature over- 
come this orbital stability, if stars are to form in the 
first place? 

The mechanism(s) responsible for generating and 
sustaining this transport of angular momentum are not 
well understood. Different options have been invoked 
by many authors over several decades, with various de- 
grees of success. It is known, however, that the molec- 
ular viscosity of accretion discs is far too low to explain 



the inferred accretion rates of these objects (Pringlc 



1981 see also Frank et al. ||1992 ). Specifically, viscous 
diffusion is able to propagate disturbances in timescales 
~ I 2 /v, where I is the distance the disturbance trav- 
els as a result of the action of the kinematic viscosity 



v. Taking I 
timescale is ' 



10 5 



■> 10 10 cm and v 
3 x 10 7 years (Balbus 



this 



2003), far too 



long to explain the variability observed in some accret- 
ing systems. Gravitational instabilities (e.g. Adams, 



Ruden & Shu 1989 Durisen et al. 2007) may be a 



viable transport process, but this mechanism may be 
effective mainly in the outer regions of massive discs 
or at early stages of the star formation process, when 
the mass of the disc is a relatively large fraction of that 
of the central object. Other proposed mechanisms are 
not general enough (i.e. the presence of a companion 
star) or do not work at all (i.e. convection). In recent 
times it has been realised that the answer may lie in the 
complex interaction of gas and magnetic fields present 
in the disc (iBalbus & Hawley 1119911 119981). Magnetic 



activity, in other words, may be the critical ingredient 
in the recipe nature follows to assemble new stars. 

From the above discussion, it is clear that angular 
momentum lies at the core of accretion-disc dynam- 
ics, and to understand how it is transported it is cru- 
cial to pay close attention to the physical conditions of 
the fluid. Most analytical and numerical studies con- 
ducted so far have adopted a number of simplifications 
to treat the fluid, and the interactions of matter with 
the magnetic field, in accreting systems. In the cold, 
dense environments typical of protostellar discsQ it is 

1 Note that the terms 'protostellar' and 'protoplanetary' discs are 
used interchangeably in the text. 



particularly important to adequately treat the slippage 
between the magnetic field and the gas (the magnetic 
diffusivity) that results from the weak ionisation of the 
disc matter. 

This paper is organised as follows. Section [2] briefly 
discusses the key ways in which magnetic fields may im- 
pact on disc physics, paying particular attention to the 
concept and types of field-matter diffusivity. Section 
[3] discusses the 'angular momentum transport problem' 
in accretion discs and introduces the most promising 
mechanisms to transfer it away, enabling accretion to 
proceed: the action of the magnetorotational insta- 



bility (hereafter MRI; |Balbus fe Hawley ||1991| |1998[ ) 
and the acceleration of outflows from the disc surfaces 
(|Konigl fe Pudritz ||2000} |Pudritz et a"l~||2007l |KoliigT 
fe Salmeron |2011 ). These mechanisms are described in 



more detail in Sections [4] and [5j respectively; and Sec- 
tion [6] discusses the possibility of their joint operation. 
Section [7] briefly summarises the main points discussed 
in the paper and concludes. 



2 Magnetic activity in protostellar discs 

Protostellar discs are known to be magnetised. Evi- 
dence of this comes from different sources. On the one 
hand, there is strong evidence for an enhanced magnetic 
activity in young stellar objects (e.g. see the reviews by 
|Feigelson et al. |2007| [Glassgold et aT~||2000[ ). As these 
authors point out, it is likely that the strong magnetic 
field near the stellar surface, extends out to the circum- 
stellar disc. On the other hand, the remnant magneti- 
sation of primitive meteorites suggests that magnetic 
fields were important in our own solar nebula as well 
( |Levy fe SonetF|[l978| . 

This property of discs is important for the present 
topic because magnetic fields are thought to play key 
roles in the dynamics and evolution of protostellar discs. 
As mentioned above, they redistribute angular mo- 
mentum via magnetohydrodynamic (MHD) turbulence 



(Stone et al. 2000), as well as by magnetocentrifugal 



winds accelerated from the disc surfaces and also pos- 
sibly from the star/disc magnetosphere (e.g. see the re- 
views byfKonigl fc Salmeron ||20lT| Konigl & Pudritz 
|2000| |Pudritz et al. |2007[ ). Additionally, magnetically- 
driven mixing strongly influences the chemistry of the 
disc (e.g. |Semenov et~aT~||2006[ |Ilgner fc Nelson ||2006[ ) 
and critically affects the dynamics, aggregation and 



overall evolution of dust particles (Turner et al. 2006 



Ciesla 2007) - the assembly blocks of planetesimals 



according to the 'core accretion' model of planet for- 
mation (Pollack et al. 1996). Furthermore, magnetic 



fields can also modify the response of the disc to grav- 
itational perturbations introduced by forming planets, 
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in turn influencing the rate and even the direction of 



planetary migration through the disc ( 


Terquem 2003 


Fromang et al. |2005||Muto et al. |2008 


Johnson et al. | 


2006). Finally, jets accelerated via magnetic stresses 



may be the sites of chondrule formation (primitive, 
thermally-processed pieces of rock found in meteorite 
samples; see e.g. |Uesugi et al. |2005| [Desch et al. |2010| ; 
and magnetically-driven activity near the disc surface 



can produce a hot, tenuous corona, (e.g. Fleming & 
Stone |2003 ) and influence the observational signatures 
of these objects. 

In protoplanetary discs, however, the magnetic dif- 
fusivity can be severe enough to limit - or even com- 
pletely inhibit - the above-mentioned processes. The 
particular role magnetic fields are able to play in these 
environments is, therefore, largely determined by the 
degree of coupling between the field and the neutral 
gas. This important topic is discussed next. 

2.1 Magnetic diffusivity 

The electrical conductivity of a fluid is determined by 
the abundance and drift motion of the charged species 
embedded in it. These include, in general, ions, elec- 
trons and charged dust grains. The abundance of these 
charged species - measured by the ionisation fraction - 
results, in turn, from the balance of ionisation and re- 
combination processes acting in the system under con- 
sideration. In protostellar discs, in particular, thermal 
ionisation is ineffective outside the central ~ 0.1 AU 
around the protostar, as the temperature drops below 
~ 1000 K. In the bulk of the disc, therefore, the most 
active ionising agents are stellar X-ray and UV radia- 



tion, interstellar cosmic rays (see e.g. 


Hayashi 


1981 


Semenov et al. 2004 


Glassgold et al. |2005 


), and - to a 



much lesser extent - the decay of radioactive elements 
present within the disc (primarily 40 K; see Umebayashi 



fc Nakano |1981[|1990| ). The attenuation column of X- 
rays and cosmic rays is 10 and 100 g cm -2 , respectively. 
In contrast, according to the minimum- mass solar neb- 



ula model (Hayashi 1981), the surface density of the 
disc at 1 AU is ~ 1700 g cm -2 . Therefore, the result- 
ing fractional ionisation in the disc interior, specially 
for radii < 10 AU, is unlikely to be high enough to pro- 
vide good coupling between the charged and neutral 
components of the fluid over the entire cross-section of 
the disc. 

The degree of coupling between the charged species 
and the neutral gas in diffusive media is typically mea- 
sured by the Hall parameter f3j, the ratio of the gy- 
rofrequency to the collision frequency of species j with 
the neutrals. Here, we specialise to ions and elec- 
trons, and denote them by subscripts i and e, respec- 



tive!^] This parameter measures the relative impor- 
tance of the Lorentz and drag forces on the motion 
of the charged species. The magnetic diffusivity is a 



tensor (Cowling 1976 Norman & Heyvaerts 1985 
Nakano fc Umebayashi ||1986 Wardle ||1999 ) whenever 
the gyrofrequency of some charge carriers is larger than 
their frequency of collisional exchange of momentum 
with the neutrals (or \/3\ > 1 for at least some of the 
charged species). This is expected, and just reflects the 
anisotropy that the magnetic field is able to impart on 
the drift of the charges when at least some of these are 
well coupled to it. In such environments, the diffusivity 
tensor components are the Pedersen (?yp), Hall (ija) and 
Ohm (?7o) terms. We also often use the ambipolar dif- 
fusivity term, given by t]a = r/p — r/Q . When the charge 
carriers are ions and electrons, it can be shown that the 
following relations apply (Wardle 2007): rjn = PeVO 
and r]A = PiPeVo- Note also that 770, r]A and rjp arc 
always positive (as 770 is independent of B, whereas t]a 
depends on B 2 ), but rfn scales linearly with B and can, 
therefore, be either positive or negative depending on 
the field polarity (e.g. the sign of B z ). 

The relative drifts of the charged species with respect 
to the bulk of the neutral gas delineate three different 
conductivity regimes, as detailed below. 

• Ambipolar diffusion limit. In this regime, which is 
important at relatively low densities, most charged 
species are primarily tied to the magnetic field 
through electromagnetic stresses (e.g. \p c \ 3> \P\\ 3> 
1, or r/A ^> |?7h| 3* Vo)- This implies that the mag- 
netic field is effectively frozen into the ionised com- 
ponent of the fluid and drifts with it through the 
neutrals. This diffusion mechanism is thought to 
be an important contributor to the loss of magnetic 
support in molecular clouds, facilitating the onset of 
cloud collapse. In protostellar discs, this regime is 
expected to be dominant at radial distances beyond 
~ 10 AU and near the surface closer in. 

• Resistive ( Ohmic ) limit, in which the ionised species 
are mainly tied to the neutrals via collisions (1 3> 

IAjI > or > |»7hI > Va)- Tnis regime is 
predominant in high-density environments, typically 



2 For charged species j, endowed with mass rrij and charge \Zj\e, 
jij = (\Zj\eB/m,jc)(l/fjp). The first parenthesis is the gy- 
rofrequency around a magnetic field of signed magnitude B = 
[B| sgn{B z } (so that the j3j retain their dependence on the mag- 
netic field polarity); and c is the speed of light. The second paren- 
thesis measures the collision frequency of species j with the neu- 
tral gas, of density p. The quantity -y» = (av)j/(mj +m), where 
{crv)j is the rate coefficient of momentum exchange of species j 
with the neutrals, and m = 2.33 mjj is the mean mass of the 
neutral particles in units of the hydrogen mass. 
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Fig. 1 Magnetic diffusivity regimes - Ambipolar, Hall and Ohm - in a log nu [cm -3 ] - log B [G] plane for T = 280 
K when the charged particles are ions and electrons only, in weakly-ionised gas (see also Wardle | p007). The different 
regions are delineated by the values of the ion and electron Hall parameters |/3i| and \0 e \, respectively. The /3j measure the 
ratio of the gyrofrequency to the collision frequency of charged species j with the neutrals and are, therefore, a measure of 
the relative importance of the Lorentz and drag forces on the motion of the charged species (see text for further details). 
The ions (electrons) decouple from the magnetic field when |/3i| (|/3 G |) drops below unity. This occurs to the right of the 
appropriate inclined line in the figure. 



close to the midplane in the innermost regions of pro- 
tostellar discs (R ~ 0.1 — 1 AU), just outside the re- 
gion in the vicinity of the protostar where thermal 
ionisation is effective. In this limit, the collision fre- 
quency of all the charged species with the neutrals 
is high enough to effectively suppress the drift of the 
former through the bulk of the gas. 
Hall limit, characterised by a varying degree of mag- 
netic coupling amongst charged species. This regime 
is important at intermediate densities, in-between 
those at which the ambipolar diffusion and resistive 
limits are dominant. In this regime, some charged 
species (typically electrons) are still tied to the mag- 
netic field whereas more massive particles (such as 
ions and charged dust grains) have already decou- 
pled to it, and are collisionally tied to the neutral gas 
(|/3 C | > 1 > IAD- In this limit, therefore, |t/h| > Va 
and rjo- As discussed above, when this regime domi- 
nates, the magnetic response of the disc is no longer 
invariant under a global reversal of the magnetic field 
polarity (e.g. |Wardle k Ng |1999||K5nigl et al. |2010 



hereafter KSW10). The Hall regime is expected to 
dominate over fluid conditions satisfied over an ex- 



tended range of radii in protostellar discs (e.g. Sano 
k Stone ||2002a|b[ ) 



In the collisionally-dominated resistive limit, the re- 
sulting magnetic diffusivity is a scalar, the familiar 
Ohmic resistivity. This is simply because collisions oc- 
cur in all directions, therefore the impulse they commu- 
nicate to the charges is randomly oriented. In the other 
limits, however, where at least some charged species are 
well tied to the magnetic field, the diffusivity is a tensor. 

Previous studies on the magnetic activity of discs 



(e.g. Wardle 


1999 Sano k Stone |2002a|b Salmeron k 


IWardle |2003 


2005| Wardle & Salmeron |2011 hereafter 



WS11) have highlighted the importance of incorporat- 
ing in these studies all the three field-matter diffusion 
mechanisms described above, as their relative impor- 
tance is a strong function of location within the disc. 
This is illustrated in Fig. [T] which shows the regions 
where the ambipolar, Hall and Ohmic diffusivity terms 
dominate in a Log B - Log «h plane (e.g. Wardle 



2007 1 , assuming that the charged particles are ions and 
electrons, in a weakly-ionised plasma. As expected, am- 
bipolar diffusion is dominant at low densities and strong 
fields whereas the opposite is true for the Ohmic regime. 
The large, intermediate region of parameter space be- 
tween these two limits is dominated by the Hall diffu- 
sivity. For example, for r = 1 AU the number density 
at the disc midplane, according to the minimum-mass 
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solar nebula model ( Hayashi 1 1981 1, is ~ 6 x 10 14 cm -3 . 
This implies that, under the adopted assumptions, the 
gas at this location would be in the Ohmic diffusivity 
regime for field strengths < 2.5 G and in the Hall limit 
for stronger fields. 

More generally, the location of the curves delineat- 
ing the regions where each of the diffusivity regimes 
dominates is a function of the mass of the charged par- 
ticles and the fractional ionisation of the fluid. If the 
charged particles are more massive, the Ambipolar-Hall 
and Hall-Ohm transitions will occur at stronger fields 
for a given value of the neutral density. Similarly, if the 
ionisation fraction increases (e.g. the neutral density is 
a lower ratio of the total density of the fluid), the above 
transitions will take place at comparatively lower mag- 
netic field strengths. A particularly important scenario 
is when dust particles are well mixed with the gas. 

The degree of field-matter coupling is strongly de- 
pendent on the abundance and size distribution of dust 
particles mixed with the gas. Dust grains tend to dra- 
matically reduce the fractional ionisation of the fluid, 
as charged particles can stick to them and recombine 
on their surfaces. Additionally, grains can also become 



an important species in high-density regions (e.g. Ume- 



bayashi fe Nakano 1990 Nishi et al. 1991). As they 



have relatively large cross sections, grains typically be- 
come decoupled from the magnetic field at lower den- 
sities than those for which other -smaller- species do. 
As a result, the conductivity of the fluid diminishes 
when dust particles (particularly when these are small) 
are well mixed with the gas. This is most relevant at 
early stages of the accretion process, or when turbulent 
motions prevent the grains from settling towards the 
disc midplane. In fact, recent diffusivity calculations 



for a minimum- mass solar nebula disc at 1 AU ( Wardle 



2007| WS11) have shown that the magnetic coupling 



may be adequate over the entire disc cross-section, pro- 
vided that dust grains are settled out of the gas phase. 

A word of caution applies to the above statements. 
The turbulent stirring of dust grains is, evidently, de- 
pendent on the disc being able to sustain turbulence 
when dust grains are well mixed with the gas. How- 
ever, as disc turbulence is thought to be a magnetically- 
driven process (see Section [5]), the reduction in the de- 
gree of field-matter coupling mediated by dust grains 
may, in turn, weaken the turbulent motions, or even 
prevent their development in dusty plasmas. Disc dy- 
namics is then the result of a complex interplay: When 
dust grains are well mixed with the gas, the disc may 
not be susceptible to turbulence, and the grains may 
efficiently settle to the midplane. On the other hand, 
once the grains are effectively removed from the gas, 
turbulence may be viable again, and the turbulent mo- 
tions may stir the grains back into the gas phase. The 



resulting structure and evolution of discs, incorporat- 
ing this important feedback effect, have not yet been 
calculated self-consistently. 



3 Angular momentum transport 

Most current models of accretion discs are viscous mod- 
els that adopt, in some form, the a-prescription of 
Shakura & Syunyaev (1973). In this formulation, the 



radial- azimuthal component of the stress tensor is as- 
sumed to scale with the gas pressure, or w r( p = aP. 
The associated turbulent viscosity is parameterised as 
v t = ac s h^, where c s is the isothermal sound speed 
and /it is the tidal density scale-height of the disc 
(e.g. resulting from the vertical gravitational compres- 
sion) . Despite its usefulness in providing an operational 
framework to support the modelling of accretion discs, 
this methodology does not offer any explanation for the 
physical nature of the accretion torque, whose origin re- 
mains unspecified. In fact, in this formulation, all the 
unknowns and uncertainties relating to the accretion 
stress are just effectively lumped into the convenient 
parameter a. 

In general, one of the most natural ways of generat- 
ing turbulence is via hydrodynamical instabilities. Lab- 
oratory shear flows, for example, readily break into tur- 
bulence at sufficiently high Reynold^ numbers. How- 
ever, Keplerian disc^] satisfy the Rayleigh's hydro- 
dynamical stability criterion (angular momentum in- 
creases with radius), so the generation and sustaining 
of hydrodynamic turbulence in these systems is, at best, 
unproven (e.g. see review by |Balbus fe Hawley |p98| . 
Convective turbulence has also been considered as an 
option ( Lin fc Papaloizou~||1980 Lin et al7] [l993), but 
further studies appear to indicate that this mechanism 
transports angular momentum towards the central ob- 



ject ( 


Cabot & Pollack 


1992 I 


lyu & Goodman 


1992 


Cabot 


1996 


Stone & Balbus 


1996 


). More generally, 



the sign of the radial flux generated by convective tur- 
bulence may depend on the ratio of the epicyclic fre- 



quency and the wave frequency of the fluid (Balbus 
2003): Angular momentum is likely to be transported 
inwards by long-period, incompressible perturbations; 



3 The Reynolds number R c = vl/u, where v is the flow velocity, 
v is the kinematic viscosity and I is an adopted length scale of 
the flow. i?, c then represents the ratio of the inertial to viscous 
forces in the fluid. 

4 In a Keplerian disc the orbital velocity is given by the relation 
% = y/ GM/r, where G is the gravitational constant, M is the 
mass of the central object and r is the cylindrical radius. The 
ratio J!k = vj^/r is the Keplerian angular frequency. 
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and outwards by high-frequency, compressible distur- 
bances. Still, other authors have pointed out that hy- 
drodynamic waves originated by the gravitational field 
of a companion star could possibly transport angular 



momentum 


Vishniac & Diamond 


1989 


Rozyczka & 


Spruit ||1993 


. This mechanism, however, makes use of 



an external source to excite, and maintain, such density 
waves; therefore it can not explain accretion in stars 
that do not belong to multiple systems. Finally, al- 



though gravitational instabilities (e.g. Adams, Ruden 



& Shu 1989 Durisen et al. 2007) may play a role, 



these may be effective mostly during early stages of 
evolution of these objects, and for particularly massive 
discs. 

It is increasingly clear that the two most generally 
applicable angular momentum transport processes are 
mediated by the action of magnetic fields and their in- 
teraction with disc matter. They are the following. 



1. 



Vertical transport via winds accelerated centrifugally 
from the disc surfaces ('disc winds'; e.g. see the re- 
views by|K6nigl fc Pudritz |2000||Pudritz et al. |2007| 
Konigl fc Salmeron ||2011 1 . This notion is strength- 



ened by the ubiquity of these outflows in star forming 
regions. 

Radial transport through turbulent viscosity induced 



by the magnetorotational instability (MRI; Balbus 
& Hawley 1119911 119981). The MRI, essentially, taps 



into the free energy contributed by the differential 
rotation of the disc, and converts it into turbulent 
motions that transfer angular momentum radially 
outwards. 

These processes are, in turn, likely to strongly affect 
the structure, dynamics, evolutionary path and planet- 
forming capabilities of their host discs. In the next two 
sections I describe in more detail the operation of these 
two forms of angular momentum transport in protostel- 
lar discs. I then briefly consider the possibility for their 
operation at the same radius, but at different heights 
above the disc midplane. 



4 Vertical transport via magnetocentrifugal 
outflows 

Fast, collimated winds are commonly observed in ac- 
creting astrophysical systems. It is likely that these 
outflows transport angular momentum away from the 
disc that surrounds the central object and thus, play an 
important role in regulating accretion (e.g. see reviews 
by|Konigl fc Pudritz |2000"t |Pudritz et al. |2007[ |K5nigl| 
fc Salmeron ||2011* ). The apparent correlation between 
outflow and accretion signatures in these systems (e.g. 



Cabrit et aT~|[l990l |Cabrit fc Andre |[199T] |Hartigan 
et al. ||1995 1 , lends support to this interpretation. The 
details of this inflow-outflow mechanism are far from 
being well understood, but it is thought that this pro- 



cess may be mediated by magnetic stresses Blandford 



fc Payne 11982} |Wardle fc Konigl |1993| (hereafter BP82 
and WK93, respectively) ; |Li ||1996" " 



Fig. [2] shows a typical disc-wind solution, plotted as 
a function of z/h^- This solution has been computed 
using the procedure detailed in KSW10 and Salmeron 
et al. (2011 hereafter SKW11); in turn a generalisation 
of the formulation by WK93. The model parameters 
are: 

1. a = «ao/c s , the midplane (subscript '0') ratio of the 
Alfven speed to the isothermal sound speed. This 
parameter measures the magnetic field strength; 

2. The midplane ratios of the diffusivity tensor com- 
ponents: (?7p/?7_l)o [or, equivalently, (?7h/?7jJo], an d 
( r l-i-/ T lo)o, where n± — \Jn\ + n\ is the total diffu- 
sivity perpendicular to the magnetic field; 

3. Ao = f ao/^J-O^k, the midplane value of the Elsasser 
number, which measures the degree of coupling be- 
tween the neutral fluid and the magnetic field (with 
A > 1 and -C 1 corresponding to strong and weak 
field-matter coupling, respectively); 

4. e = —v r o/c s , the normalised inward radial speed of 
the fluid at the midplane. This parameter is evalu- 
ated by imposing the Alfven critical-point constraint 
on the wind solution (see SKW11); 

5. c s /vk — hi/r, the ratio of the disc tidal density 
scale-height to the radius (e.g. the disc geometric 
thickness); and 

6. eg = — «Bro/ c s = —cE^q/ c s B z , the normalised (and 
constant with z) azimuthal component of the elec- 
tric fielcj^J This parameter measures the radial drift 
speed of the poloidal magnetic field lines (WK93). 

We will also use the expression T = /^i/Slx, the 
ratio of the Keplerian rotation time to the neutral- 
ion momentum exchange time (e.g. WK93; see also 
footnote [2]). It can be shown that in the ambipolar 
regime A = T, whereas in the Hall limit A = T|/3;| 
(e.g. KSW10). 

The left panel of Fig. [2] displays the gas density 
(normalised by its value at the midplane) as well as 
the radial and azimuthal components of the magnetic 
field (normalised by the vertical field, which under the 
adopted thin-disc approximation is vertically constant), 
as a function of height above the midplane, expressed 

5 Setting eb = effectively fixes the value of B r at the disc sur- 
face. In a global formulation, B r would be determined by the 



radial distribution of B z (e.g. |Ogilvie fc Livio |200l) |Krasnopol-| 
|sky fc Konigl |2002[ l. 
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Fig. 2 Vertical structure of a strongly-magnetised, wind-driving disc solution, as a function of height above the disc 
midplane (in units of /it, the tidal density scale- height), under the assumption of pure ambipolar diffusivity. The model 
parameters are: A = 10, ao = 0.95, ?7po£V c s = 0.09, t]hoQ/Cs = 0, e = 0.7, ep, = 0, and vk/c b = 20 (see text for details). For 
simplicity, in this illustrative solution the diffusivity components are assumed to scale with height in such a way that the 
Elsasser number (A) remains constant with z. In a more general formulation, these parameters would vary with height as 
dictated by the ionisation balance in the disc. All curves terminate at the sonic point. Left: Gas density (p/po, normalised 
by the midplane value) and magnetic field components (radial and azimuthal, B r /B z and \B$\/ B z , respectively; normalised 
by the vertical field B z , which is constant with height). In the shaded region, the gas thermal energy dominates over 
the magnetic energy and the magnetic field lines are bent and sheared. Right: Velocity components, with respect to the 
Keplerian velocity vk and normalised by the isothermal sound speed c s . Within the disc (to the left of the red vertical line), 
the radial velocity is negative, the azimuthal velocity is sub-Keplerian and the vertical velocity is small. On the contrary, 
in the outflow region above this layer all three velocity components [(v — vk^/cb] are positive and increasing with z. 
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Fig. 3 Vertical structure of a strongly-magnetised wind-driving disc (see WK93 and also Konigl fc Salmeron |2011 1. Three 
distinct layers can be identified (see text): (1) a quasi-hydrostatic layer, next to the midplane (z < Zh, where Zh s ft is 
the magnetically-reduced density scale-height), which contains the bulk of the matter; (2) a transition layer (zh < z < Zb, 
where subscript 'b' denotes the launching height or base of the wind), where the magnetic energy becomes dominant and 
the inflow gradually diminishes; and (3) an outflow region (z > z^), which constitutes the base of the large-scale wind. z s 
is the height of the sonic point and vb is the field-line velocity. 
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in units of the tidal scale- height. The right panel shows 
corresponding runs for all velocity components, mea- 
sured with respect to the Keplerian velocity i>k and 
normalised by the isothermal sound speed c s . This so- 
lution illustrates the main features of winds accelerated 
centrifugally from strongly-magnetised discs (WK93). 
Three distinct layers can typically be identified. These 
layers are presented schematically in Fig. [3j and de- 
scribed below. 

• The quasi-hydrostatic layer, adjacent to the disc mid- 
plane (z < Zh = h, where h is the magnetically- 
reduced density scale-heightQ) , is matter dominated 
(vb4> < v <f> < vk, where vb is the field-line velocity). 
In this layer, the topology of the magnetic field is 
consistent with the field lines being radially bent and 
azimuth-ally sheared. The field removes angular mo- 
mentum from the g result of the ion-neutral 
drag, and matter moves radially inwards. 

• The transition region, above the quasi-hydrostatic 
layer (z^ < z < Zb, where subscript 'b' denotes the 
base of the wind), where magnetic energy comes to 
dominate over the thermal energy of the fluid as a re- 
sult of the strong drop in density away from the mid- 
plane. The magnetic field lines are locally straight in 
this zone (note that they arc still inclined outwards). 
Here (as in the quasi-hydrostatic layer at lower z), the 
fluid azimuthal velocity is sub-Keplerian, as magnetic 
tension helps support the fluid against the gravita- 
tional pull of the central protostar. 

• As the angular velocity of the field lines is constant, 
their azimuthal velocity increases along the line, until 
they eventually overtake the fluid in which they are 
embedded. This location marks the beginning of the 
uppermost wind region (z > Zb), which constitutes 
the base of the outflow. In this layer the field trans- 
fers angular momentum back to the matter, acceler- 
ating it centrifugally; the radial and vertical velocity 
components are positive, and the azimuthal velocity 
is super-Keplerian (ub^ > v<f, > vk)- 

In the specific solution shown in Fig. [2j the quasi- 
hydrostatic layer, which encompasses the section of the 
disc associated with strong gradients in the fluid density 
and magnetic field components, extends from z/hr = 
to z//it ~ 0.4 (the region shaded in red in the left 
panel of the figure). The transition region lies directly 
above it, and extends up to the height where the az- 
imuthal velocity becomes Keplerian (up to the vertical 

6 Magnctic compression dominates over the gravitational tidal 
squeezing of the disc in all viable wind-driving disc solutions 
of interest here (WK93, KSW10), so that the magnetically- 
compressed density scale- height (h) satisfies h/hi < 1. 



red line in the right panel), so it constitutes the re- 
gion where 0.4 < z/hi < 1.8. Finally, the wind region 
satisfies z/hi > 1.8. Note that in the first two layers 
(i.e. within the disc) the radial velocity is negative, the 
azimuthal velocity is sub-Keplerian and the vertical ve- 
locity is small. In the wind layer, on the other hand, all 
these velocities [(v — vk4>)/ c s] are positive and increas- 
ing with z. 

Solutions such as the one shown in Fig. [2] are found 
for relatively strong magnetic fields, such that oq < 1. 
Under these fluid conditions, the MRI is expected to 
be suppressed because the wavelength of the critical 
MRI-unstable mode (longward of which the instabil- 



ity operates; see Balbus & Hawley 1991) exceeds the 
magnetically-reduced density scale- height h (WK93). 
For the particular solution shown in the figure, the mag- 
netic field is strongly coupled to matter even at the disc 
midplane (Ao 3> 1), and ambipolar diffusion dominates 
over the entire cross-section of the disc. For simplicity, 
the diffusivity components are assumed to scale with 
height in such a way that the local Elsasser number 
(A) remains constant with z. Under these approxima- 
tions, it is self-consistent to find relatively high values of 
the inward velocity of the fluid at the disc midplane (of 
the order of the isothermal sound speed, such that the 
parameter e is in the range of 0.3 - 1). A more realis- 
tic treatment of the fluid conditions would incorporate 
the vertical stratification of the diffusivity, with differ- 
ent regimes dominating at different vertical locations 



(see section 2.1), and possibly a magnetically weakly- 
coupled region in the disc interior. In the latter case, 
in particular, more moderate values of the inward flow 
speed at the midplane are expected ( Li ||1996| |Wardle"| 
1997| see also Fig. [5] below) . 



The next issue to consider is that in order to compute 
self-consistent, coupled disc-wind solutions, it is essen- 
tial to solve simultaneously for the complex interactions 
between the gas and the magnetic field within the disc - 
where the launching process occurs - as well as to follow 
the acceleration of the wind past the critical surfaces 
of the flow. Physical quantities of critical importance, 
such as the gas density and the ionisation fraction, are 
expected to vary by many orders of magnitude between 
the disc midplane and its surface. For example, the 
gas density typically drops by ~ 6 orders of magnitude 
across the half-thickness of the disc. As a result, global 
simulations of these winds have adopted a number of 
simplifications to treat the underlying disc. In some 
(including the pioneering work of BP82), the disc is 
treated as a boundary condition, outside of the com- 
putational domain. Other authors have attempted to 
incorporate the disc, either via assuming perfect field- 
matter coupling (i.e. negligible diffusivity) or through 
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the prescription of a simplified finite resistivity term 



Ferreira 

| |2000| iZanni < t al. 



(e.g. see the review by Pudritz et al. 2007 and also 



2007) 



Pelletier ||1995| |Li ||1996| |Casse fc Ferreira| 
This approach is appropri- 



ate to follow the propagation of the outflow to large 
distances for an essentially prescribed wind mass flux, 
but it can not be used to determine, self-consistently, 
whether the wind is launched in the first place. 

A complementary approach, first developed by 
WK93 for an ambipolar-diffusion-dominated disc, is 
to resolve explicitly the vertical stratification of the 
fluid variables, ionisation fraction and magnetic diffu- 
sivity, as well as the resulting field-matter coupling, 
for a radially-localised region of the disc. This formu- 
lation systematically reduces the nonlinear governing 
equations that control the evolution of the fluid into 
a set of ordinary differential equations (ODE) in z. 
This system describes, for a particular radial location, 
the vertical structure of the disc - and the base of the 
outflow - from the midplane up to the critical (sonic) 
surface of the flow (see also KSW10 for a generalisa- 
tion of this procedure to a tensor diffusivity) . The 
system of ODE is numerically integrated by applying 
boundary conditions at the midplane and at the sonic 
point (see WK93 and SKW11 for details of the numer- 
ical procedure). The ionisation balance is calculated 
by evaluating the action of ionisation processes (driven 
by cosmic rays, X-rays emitted by the central object 
and radioactive decay) and recombinations occurring, 
in general, both in the gas phase and on grain surfaces 
(see Section 2.1 1. The X-ray ionisation rate is taken 
from the Monte Carlo calculations of Igca & Glassgold 
[] (1999), whereas that of cosmic rays is calculated by 
attenuating the standard ("canonical") rate in the in- 
terstellar medium as the cosmic rays penetrate the disc 
[the canonical rate is taken to be 10 17 (s _1 ) per hydro- 
gen atom; see e.g. Umebayashi & Nakano 1981 . The 



field-matter diffusivity is treated as a tensor, incorpo- 
rating the Hall, Ohmic and Ambipolar terms (Section 
2.1). This approach eliminates the need to keep sepa- 



rate equations of motion for each fluid component and 
greatly simplifies the calculations, particularly when 
dust grains are present. 

The methodology just described is appropriate to 
model the launch of the wind, but it cannot be used 
to follow the propagation of the outflow far from the 
source, where the thin-disc approximation breaks down. 
It is, however, clearly necessary to ensure that the ob- 
tained 'disc solution' continues to accelerate past the 
sonic point. To ensure that this is the case, we eval- 
uate the parameters of the outflow at the disc sur- 
face and match this local solution to a global wind 
model (we use for this purpose the BP82 self-similar 



wind solutions). In order to facilitate the matching, 
we constructed a library of these self-similar 'wind so- 
lutions' for a wide range of their model parameters, 
which are: the normalised mass-to-magnetic flux ra- 
tio (k), the normalised total specific angular momen- 
tum carried away by the wind (A, which incorporates 
contributions from the matter and from the magnetic 
field), and the inclination of the field lines at the disc 
surface = B r ^/B z ). The combinations of these pa- 
rameters corresponding to viable global wind solutions 
are shown graphically in Fig. [4] The complete set of 
these results is also available in the electronic version 
of SKW11. The described overall methodology is use- 
ful to explore the viability of these winds as a function 
of radius from the central object, as well as to obtain 
a realistic estimate of the wind mass loss rate, which 
could be used to constrain theoretical models against 
observational results. 

A representative local solution, matched to the BP82 
global, self-similar models, is shown in Fig. [5] for r = 1 
AU in a disc surrounding a solar-mass protostar. The 
surface density is £ = 600 g cm~ 2 and the remaining 
free model parameters are: ao = 0.75, «k/c s = 10, 
e = 0.1 and £b = 0. The magnetic diffusivity com- 
ponents and field-matter coupling (A) are calculated 
self-consistently by evaluating the ionisation balance 
of the fluid (see Section 2.1). On the left-hand side 



of the figure, the top panel shows the ionisation rates 
contributed by cosmic rays (curve labelled by the sub- 
script 'cr'), X-rays ('xr') and radioactive decay ('rad'); 
and the bottom panel displays the resulting diffusivity 
terms. On the right-hand side, the normalised density 
and magnetic field components in the plane of the disc 
are shown in the top panel. The normalised velocity 
components are displayed in the bottom panel. The 
global wind model parameters are: k = 2.6 x 10 -6 , 
A = 4.4 x 10 3 , and C h = B rh /B z = 1.6. The local mass 
accretion rate for this model is ~ 7 x 10~ 6 Mq yr _1 , 
which is consistent with the inferred values for the early 



(Class 0/Class I) protostellar accretion phase (e.g. Har- 
tigan et al.~||1995[ ) 



We have also constrained the parameter space oc- 
cupied by physically-viable wind-driving disc solutions 
(see WK93 for the original derivations in the ambipolar 
diffusion limit; as well as KSW10 for the generalisation 
to a tensor diffusivity and application to the Hall and 
Ohm regimes). Using the hydrostatic approximation 
(e.g. neglecting the vertical velocity component of the 
fluid), we are able to show that the solutions are re- 
quired to satisfy the following constraints. 

• Sub-Keplerian flow below the launching region. Be- 
low the launching height the azimuthal velocity of 
the gas should be sub-Keplerian, as discussed earlier 
in this section. 
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1 2 3 4 5 

Log(A) 

Fig. 4 Global, self-similar wind solutions of BP82 in parameter space. The solutions are plotted in a £{, (= B r b/B z ) - log 
A plane, where A is the normalised total specific angular momentum carried away by the wind (including the matter and 
the magnetic field contributions) and £{, measures the inclination of the field lines at the base of the outflow. The curves 
are labelled by the normalised mass-to-magnetic flux ratio (ft) and the values of k that correspond to the darker curves are 
indicated in the figure. See SKW11 for the method of calculation of these solutions and for a tabular compilation of the 
results. 
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Fig. 5 Representative radially-localised disc-wind solution at r = 1 AU, matched to the global, self-similar models of 
BP82, for a disc surrounding a solar-mass protostar. The surface density of the disc is E = 600 g cm -2 and the remaining 
free parameters are: ao = 0.75, vk/c s = 10, e — 0.1 and eb = 0. The magnetic diffusivity components and field-matter 
coupling (measured by the Elsasser number) are calculated self-consistently by evaluating the ionisation balance in the disc 
(see text). The parameters of the self-similar model are k = 2.6 x 10 -6 , A = 4.4 x 10 3 , and = B r ^/B z = 1.6. Left 
column. Top panel: Ionisation rates contributed by cosmic rays (cr), X-rays (xr) and radioactive decay (rad). Bottom panel: 
Ambipolar, Hall and Ohm diffusivity components. Right column. Top panel: Density and magnetic field components in 
the plane of the disc (radial and azimuthal). Bottom panel: Velocity components. The resulting local mass accretion rate 
of the model is 7 x 10 -6 Mq yr _1 , a value that is consistent with the accretion rates inferred for objects in the early stages 
of the protostellar accretion phase (e.g. Class 0/Class I objects; e.g. Hartigan et al. |1995 |. 
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Wind-launching criterion. The sufficient inclination 
of the field lines with respect to the rotation axis of 
the disc for centrifugal acceleration to occur (BP82). 
Location of the base of the wind. Both theoretical 
and observational arguments suggest that only the 
uppermost section of the disc should participate in 
the outflow (|Konigl fc Pudritz |[2000|. We apply 



this condition by requiring that the base of the wind 
be located above one magnetically-reduced density 
scale-height, or z^/h > 1. It can be shown that when 
this condition is violated, the azimuthal component 
of the magnetic field changes sign within the disc, in- 
dicating an attempt by the field to transfer angular 
momentum back to matter below the wind launching 
height. 

• Dissipation rate. The rate of heating by Joule (e.g. 
magnetic) dissipation at the disc midplane is re- 
stricted to be less than the rate of gravitational po- 
tential energy release at that location ( Konigl |1997 ). 

The above constraints can be shown to lead to four 
parameter sub-regimes where viable solutions exist in 
the Hall regime, and to three in the Ohm regime (I refer 
the reader to KSW10 for more details). The constraints 
for the Hall regime are summarised in Table [TJ These 
are differentiated by the values of the combinations 
/3coAo an d 2Ao (= 2To|/3jo| in this limit) in comparison 
with unity. Our numerical solutions (SKW11) confirm 
these predictions: We find that no physically-relevant 
solutions are found outside the boundaries specified by 
these constraints. 

We also explored the dependence of the solutions on 
the model parameters. We found that increasing the 
relative contribution of the Hall diffusivity results in a 
smaller magnetically-reduced density scale-height, and 
in a higher location above the disc midplane of both 
the base of the wind and the sonic surface. As a re- 
sult, the density at the sonic point, and the associated 
mass outflow rate, are reduced. Our calculations also 
show that, in all diffusivity regimes, viable solutions 
satisfy the requirement that the neutral - ion momen- 
tum exchange time be shorter than the disc orbital time 
(T > 1). This is, therefore, a fundamental constraint 
on the wind solutions considered here. 

Finally, we found that the magnetic field polarity 
modifies both the properties of the solutions and the pa- 
rameter ranges where these exist when the Hall regime 
is dynamically important. This reflects the dependence 
of the Hall diffusivity on the sign of B z . Specifically, 
our analysis shows that for two of the Hall sub-regimes 
(Cases i and Hi in Table [TJ, positive and negative po- 
larity solutions have different properties: For the same 
values of all the other parameters, both the base of the 
wind and the sonic point are located higher above the 
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Fig. 6 Normalised sonic-point density (p s /po) as a func- 
tion of the parameter To for solutions in the Hall sub- 
regime (i) (KSW10, see also Table [TJ. These results illus- 
trate the dependence of physically-viable wind-driving disc 
models on the orientation of the magnetic field when the 
Hall diffusivity dominates. The red curves correspond to 
f]H = r/H^/Cs > (or B z > 0), and the blue ones to fju < 
(B z < 0). These are labelled by the normalized midplane 
inflow velocity e. The diffusivity values listed at the top of 
the figure refer to the disc midplane. The value of |/3io| _1 
is the same (1.45) for all the solutions. The results are 
broadly consistent with the prediction of the hydrostatic 
analysis (KSW10) that there should be no viable solutions 
for fju < when To drops below 2|/3io| _1 (i-e. to the left of 
the vertical dashed line). Where solutions with both polar- 
ities exist, their properties are different, as predicted. 
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Table 1 Parameter constraints for wind-driving disc solutions in the limit where the Hall diffusivity dominates. Four 
distinct cases can be identified, depending on whether the combinations /3 e o/3io and 2Ao (= 2Toj/3io|) are larger or smaller 
than unity (KSW10). The first inequality expresses the requirement that the disc remain sub-Keplerian below the launching 
point, the second is the wind launching condition (the requirement that the magnetic field lines be sufficiently inclined with 
respect to the angular velocity vector of the disc for centrifugal acceleration to occur), the third ensures that the base of the 
wind is located above the magnetically-reduced density scale-height, and the fourth specifies that the rate of Joule heating 
at the midplane should not exceed the rate of release of gravitational potential energy at that location. 
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midplane when the magnetic field is aligned with the 
angular velocity vector of the disc (B z > 0) than when 
it is pointing in the other direction (B z < 0). The 
mass outflow rate is correspondingly lower in the for- 
mer case than in the latter. Furthermore, no viable so- 
lutions are predicted to exist in these sub-regimes when 
T < 2 1 /?io I 1 an d B z < 0. For the other two parameter 
sub-regimes (Cases ii and iv in Table [l]), the situation 
is even more drastic: Solutions are predicted to exist 
for the positive polarity only (e.g. no viable solutions 
are expected for the negative polarity). 

Our numerical solutions confirm the above expecta- 
tions. This is illustrated in Fig. [6] (for sub-regime i). 
This figure shows the normalised density at the sonic 
point {p s / po) as a function of the parameter To, for so- 
lutions with both (positive and negative) magnetic field 
polarities. In all cases, the value of lAol -1 = 1-45. The 
displayed results confirm the prediction of the hydro- 
static analysis (KSW10) that no viable solutions exist 
for % = ?7h^/c 2 < when T < 2|Aol _1 (i-e. to the left 
of the vertical dashed line). It is also evident that in the 
region of parameter space where solutions with both po- 
larities exist, these have different properties. Note that 
the mass outflow rate (measured by p s / Po) is smaller in 
the B z > case, in accordance with our predictions. 



5 Radial transport via MRI-induced 
turbulence 



The MRI (Balbus fc Hawley 



view by Balbus fc Hawley |1998 ) , transfers angular mo- 



1991 see also the re- 



mentum radially outwards via magnetic field lines that 
connect fluid elements located at different radii. The 
physical principle that underpins this mechanism is il- 
lustrated in Fig. [7j Imagine two fluid elements that 



are initially located at the same orbital radius and are 
joined by a vertical magnetic field line threading the 
disc (left panel). In this configuration the magnetic 
tension in the line is negligible. Now imagine that the 
elements are perturbed from their initial orbits, so that 
one of them (the one labelled mi in the figure) is dis- 
placed to an inner orbit, whereas the other (777,2) moves 
to an outer orbit (middle panel). Assuming that the 
field-matter coupling is sufficient for the magnetic field 
line to be dragged by the moving fluid elements, the 
line will be stretched and magnetic tension will develop. 
Note that the element in the inner orbit (mi ) overtakes 
the one in the outer orbit (777,2), as a result of the de- 
pendence of the Keplerian orbital motion with radius 
(uk = \J GM/r). The tension in the line, therefore, 
acts to reduce the angular momentum of the inner ele- 
ment, and to increase that of the outer element (right 
panel). As a result of this process, angular momentum 
is transferred from mi to m 2 , causing the former ele- 
ment to move further inwards while the latter retreats 
to an even larger radius. As the elements continue to 
separate, the tension in the magnetic field line joining 
them increases, and the process "runs away" ; pertur- 
bations, once initiated, will amplify. 

The nonlinear stages of the instability have been 



studied extensively (e.g. see the reviews by Balbus & 



Hawley 2003 1998). The result is MHD-driven tur- 



bulence where angular momentum is transported ra- 
dially outwards. This outward transport is linked to 
the fact that the turbulence is anisotropic: the radial 
and azimuthal components of the perturbations of the 
magnetic and velocity fields are significantly correlated. 
In other words, the transport of angular momentum 
follows directly from the statistical correlation in the 
Reynolds (velocity) and Maxwell (magnetic field) stress 
tensors ( |Balbus fc Hawley ||1998[ ). 
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Under ideal-MHD conditions (negligible magnetic 
diffusivity) , MRI- unstable modes grow in weakly- 
magnetised discs (meaning here that the magnetic en- 
ergy density is small in comparison to the thermal en- 
ergy density; or a = Vx/c s <C 1), provided that the an- 
gular velocity of the fluid decreases radially outwards. 
This condition is satisfied in Keplerian discs, for which 
= v K/ r fx r~ 3 / 2 - An analysis of the dispersion 
relation, linking the wavenumbers and growth rates of 
the MRI-unstable modes (e.g. Balbus fe Hawley |1991 ), 
shows that the maximum growth rate of the instability 
is of the order of SIk, and it is independent of the field 
strength. 

The upper limit on the range of magnetic field 
strengths for which the MRI is expected to grow in discs 
(under ideal-MHD conditions) reflects the fact that the 
critical wavelength (longward of which the instability 



and electrons only (e.g. dust grains have settled out of 
the gas phase). They were obtained using the proce- 



operates) is directly proportional to B ( Balbus & Haw- 



ley 1991). For a sufficiently strong field, therefore, 



the critical wavelength will exceed the thickness of the 
disc (~ /it) and the instability will be damped. It is 
also the case that the wavelength of the most-unstable 



mode is larger than the critical wavelength (jBalbus fc ?7h^k/«a 



Hawley 1991). Therefore, even if the wavelength of 



maximum growth does not fit within the disc, there 
may be slower-growing modes that are able to develop 
(provided that the disc thickness is larger than the crit- 
ical wavelength of the instability). In these conditions, 
MRI-induced turbulence is still viable, although at a 
reduced strength. Because of the general conditions 
under which it operates, the MRI is a promising mech- 
anism for the redistribution of angular momentum in 
astrophysical accreting systems. 

Ideal-MHD conditions, however, are not likely to 
be satisfied in the deep, dense interiors of protostellar 
discs. In these environments, therefore, it is essential 
to consider how the high magnetic diffusivity impacts 
on the development and properties of the instability. 
Fig. [8] illustrates this point. This figure shows the ver- 
tical structure and the linear growth rate of the fastest- 
growing MRI perturbations for a minimum-mass solar 
nebula disc at r — 1 AU, as a function of the mag- 
netic field strength and for different diffusivity limits 
(Ambipolar diffusion, Hall diffusion with B z > 0, and 
the 'Full case', where both components of the diffu- 
sivity tensor are present). In each panel, the solid 
(dashed) line denotes the radial (azimuthal) perturba- 
tion of the magnetic field, which is initially vertical. 
The field strength is indicated in the top-right corner 
of each panel, whereas the growth rate, in units of Qk, 
is shown in the lower-right corner. This notation is ob- 
served here in all plots that display the structure of the 
instability. These solutions incorporate a realistic ioni- 
sation profile, assuming that charges are carried by ions 



dure detailed in Salmeron & Wardle (2003 2005), to 



which I refer the reader for additional details. Note 
that Hall diffusion strongly modifies the structure and 
growth of the MRI-unstable modes. Specifically, when 
this diffusion mechanism is included in the calculations 
the vertical extension of the magnetically-inactive zone 
straddling the disc midplane is significantly reduced, 
and the growth rates of the perturbations increase in 
comparison to those obtained in the ambipolar diffu- 
sion limit. Note that for the field strengths shown, the 
perturbations in the Hall limit grow at about the ideal- 
MHD rate, whereas the ones in the ambipolar diffusion 
limit grow at a significantly reduced rate. 

WS11 have conducted a local linear analysis of the 
MRI, for a weak vertical field and considering axisym- 
metric perturbations with a solely vertical wavevector. 
This study demonstrated that for a given value of the 
Pedersen diffusivity (rip), the maximum growth rate 
increases with the Hall diffusivity term (??h)- In par- 
ticular, for a given r/p the growth rate vanishes when 



< —2, and it increases to the ideal-MHD 
maximum of ~ 0.7517k as 7/h — > oo. The authors also 
show that, for the minimum-mass solar nebula model 
and at r — 1 AU, Hall diffusion modifies the extent of 
the magnetically active column by one order of mag- 
nitude, either increasing it or decreasing it depending 
on whether the magnetic field is parallel or antiparallel 
to the angular velocity vector of the disc, respectively. 
This trend has not yet been reflected in existing numer- 



& Stone 



ical simulations of MRI-induced turbulence (e.g. Sano 
2002a|b| . We argue that this is the result of 



these simulations not yet probing the Hall-dominated 
region of parameter space, where these effects are ex- 
pected to be significant (WS11). 

Fig.[9]compares the growth rate of the most-unstable 
MRI modes as a function of the magnetic field strength, 
for different radii in a minimum-mass solar nebula disc 
and assuming that ions and electrons are the only 
charge carriers. Note that in all cases the growth rate 
initially increases with B, then levels at u max ~ 0.75 
(the growth rate associated with ideal-MHD perturba- 
tions in a Keplerian disc) and, finally, drops abruptly 
at the field strength for which the wavelength of the 
perturbations is of the order of the tidal scale-height of 
the disc. This generally occurs when the ratio of the 
magnetic pressure to gas pressure at the midplane is of 
order unity (or a Q < 1; Balbus fc Hawley "||1991 ). How- 
ever, calculations by Salmeron & Wardle (2003) and 



WS11 show that when 7?h < the instability may grow 
for super-equipartition (supra-thermal) fields. This is 
the result of Hall diffusion in this case extending the 
unstable wavelengths to shorter values. 
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1. Suppose that two fluid 2. If they are perturbed 



elements, in the same 
orbit, are joined by a 
field line (B) 



B 



6 rn 1 



to different orbits 



3. The tension acts to 
reduce the angular 
momentum (L) of 
m t and increase 
that of m, 



AsL OC V L 



m 1 moves inwards 




The tension in the 
line is negligible 



The field line is stretched 
and develops tension 



This increases the tension and the process "runs away" 



Fig. 7 Physical principle underpinning the action of the magnetorotational instability (MRI). The left panel shows a 
schematic diagram (in the r-z plane) of a protostellar disc threaded by a vertical magnetic field. Two fluid elements are 
located at the same orbital radius, joined by the field line. This panel corresponds to the initial configuration of the disc. 
Once perturbations are initiated (middle panel), these will amplify (right panel, showing the disc viewed from the top). See 
text for details. 



Full case Ambipolar Diffusion Hall Diffusion 




z/h T z/h T z/h T 

Fig. 8 Vertical structure and linear growth rate of the most-unstable MRI modes for a minimum-mass solar nebula disc 
at r = 1 AU, as a function of the magnetic field strength (shown in the upper-right corner of each panel) and for different 
diffusivity limits. In each panel, the solid (dashed) line denotes the radial (azimuthal) perturbation of the magnetic field, 
which is initially vertical. The growth rate (in units of the Keplerian frequency Ok) is indicated in the lower-right corner 
of each panel. The middle and right columns show solutions obtained under the Ambipolar and Hall diffusion (B z > 0) 
approximations, respectively. The left column depicts solutions incorporating both diffusivity components. These results 
show that MRI perturbations grow faster, and are active closer to the midplane, when both diffusion mechanisms are 
considered. 
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These results demonstrate that the MRI grows over a 
wide range of magnetic field strengths once dust grains 
have settled out of the gas phase, even in the inner re- 
gions of weakly-ionised discs. Note also that the maxi- 
mum field strength for which perturbations grow is in- 
versely related to the radius. This is expected, as both 
the midplane gas density po and the isothermal sound 
speed c s decrease with the radial location. As a result, 
the value of a = i3 /[(47r/9o) 1 ^ 2 c s ] associated with a 
particular field strength increases with r, causing the 
instability to be damped at progressively weaker fields 
as the radius increases. 

The magnetic activity of discs is also strongly de- 
pendent on the presence and size distribution of dust 
grains mixed with the gas. This is the result of rapid 
recombinations taking place on their surfaces, as well as 
their relatively large cross-sections, which causes them 
to decouple from the magnetic field at densities for 
which other - smaller - particles are still well tied to 
it. The vertical structure and growth of the fastest- 
growing MRI modes as a function of the presence and 
size of dust grains suspended in the gas phase is illus- 
trated in Fig. [TO] for the minimum-mass solar nebula 
model at r — 5 AU. The leftmost column of the fig- 
ure shows the magnetic field perturbations when the 
grains have settled out of the gas phase, or have aggre- 
gated to sizes significantly larger to those shown. The 
other columns display the corresponding results when 
a population of single-sized grains, of the radius (a g ) 
indicated at the top of each column and constituting 
a constant (one percent) fraction of the mass of the 
gas, remains suspended with the fluid. All the mag- 
netic diffusivity components have been incorporated in 
the calculations and B z > 0. The lowest row shows 
the maximum field strength for which the perturba- 
tions grow (for a g = 1 and 0.1 /im), or the results for 
B = 100 mG (for the no grains, or a g — 3 /im cases). 
Note the extended magnetically-inactive zone adjacent 
to the disc midplane when the grains are small (sub- 
micron sized). In this case, severe magnetic diffusion in 
the disc interior prevents the magnetic field from cou- 
pling to the gas below ~ 2.5 scale-heights. When the 
grains have aggregated to about a micron in size, how- 
ever, the disc midplane remains magnetically coupled 
for all the field strengths for which the instability grows 
(B < 100 mG, see the column for a g = 1 ^m in Fig. 10 ). 



In contrast, the MRI grows for B < 1 G when no grains 
are present at this radiu^] ( |Salmeron fc Wardle [2005 ) . 
These results hold when the magnetic field is aligned 
with the rotation axis of the disc, rendering 77H > 0. 



When the field is counter-aligned to it, however, only a 
small fraction of the disc is active even after the grains 
aggregate, unless the field is in the range of 20 - 80 G 
(WS11). 

Fig. |TT] displays the growth rate of the fastest- 
growing MRI perturbations at r — 5 AU, for different 
assumptions regarding the presence and size of dust 
grains mixed with the gas and B z > 0. The range of 
field strengths over which the MRI operates is smaller 
as the grain size diminishes. This is expected, given the 
reduction in magnetic coupling for a particular B as 
the dust particles are smaller. Note, however, that the 
stabilising effect associated with the presence of small 
grains is not strong enough to completely suppress the 
perturbations, even when small particles (0.1 /im in 
size) are present. In this case unstable modes are still 
able to grow over a finite section of the disc away 
from the midplane and for a restricted range of field 
strengths. 

All the results discussed so far in this section were 
obtained assuming that the surface density profile of the 
disc is that of the minimum-mass solar nebula model, 
£(r) = 1700 rj^j 2 g cm~ 2 , where r^u is the radial dis- 
tance from the central object measured in astronomical 
units. This yields a surface density of ~ 150 g cm -2 
at 5 AU and ~ 50 g cm" 2 at 10 AU. However, the 
surface density profile in real discs may well be dif- 
ferent from these inferred values. In fact, recent ob- 



servations by Kitamura et al. (2002) and Andrews & 



Williams (20071 seem to imply that the surface den- 



sity declines more gradually with radius than predicted 
by this model. Moreover, the actual surface density in 
the disc inner regions may be smaller than that of a 



minimum- mass solar nebula disc ( Andrews & Williams 



2007). This would facilitate a deeper penetration of 



the ionising sources into the disc, and would result in 
an increased ionisation fraction closer to the midplane. 
This effect would, in turn, modify the properties of 
MRI-unstable modes in this region and, more gener- 
ally, the presence and configuration of the magnetically 
dead zone in the disc interior. On the other hand, pro- 
tostellar discs may have surface densities roughly up to 
~ 50 times the values implied by the minimum-mass 
solar nebula model without becoming gravitationally 
unstable. An increased surface density will result in a 
more extended section of the disc interior being shielded 
from the ionising effect of X-rays and cosmic rays, but 
the magnetically active regions closer to the surface are 
not likely to be substantially modified. 



7 The MRI grows for B < 8 G (at 1 AU) and B < 250 mG (at 10 
AU) when dust grains are absent (ISalmeron & Wardle 12005b. 
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Fig. 10 As per Fig.[8]as a function of the presence and size of dust grains suspended in the gas phase, for a minimum-mass 
solar nebula disc at r = 5 AU. All diffusivity components have been incorporated in the ionisation balance, and B z > 0. The 
leftmost column shows the magnetic field perturbations when the grains have settled out of the gas phase. The remaining 
ones display (from left to right) the corresponding results when a population of single-sized grains of radius a g = 3, 1 and 
0.1 /im, respectively, remains suspended with the fluid. The grains constitute a constant fraction (one percent) of the total 
mass of the gas. The lowest row shows the maximum field strength for which the perturbations grow (for a s — 1 and 
0.1 lira), or the results for B = 100 mG (for the no grains and the a s = 3 /im cases). Note that when the grains are small 
(a g = 0.1 /zm) the layer below ~ 2.5 scale-heights is magnetically inactive. When B z > 0, however, this effect is restricted 
to a finite range of grain sizes. Once the grains have aggregated to ~ 1 (im or more, the midplane is again unstable to the 
MRI (see also |Salmeron fc Wardle |2008| and WS11). 
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Fig. 9 Linear growth of the most-unstable MRI modes for 
a minimum-mass solar nebula model at r = 1, 5 and 10 AU 
as a function of the strength of the magnetic field. For each 
radius, the maximum field strength for which the instability 
grows is the one for which the ratio of the Alfven speed to 
the isothermal sound speed at the midplane is of order unity 
(ao ~ 1)- For this B, the critical wavelength (longward of 
which the MRI operates) is of the order of the disc thickness 
(~ /it) and the instability is damped. 
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Fig. 11 As per Fig. [9] for r = 5 AU and different as- 
sumptions regarding the presence and size of dust particles 
suspended in the gas. The range of field strengths over 
which the MRI operates is more restricted as the grains are 
smaller, given the associated drop in the field-matter cou- 
pling for a particular B as the grain size goes down. 



6 Combined radial and vertical angular 
momentum transport 

The redistribution of angular momentum in real proto- 
planetary discs is likely to involve both the radial and 
vertical transport mechanisms discussed in the previous 
two Sections. Although there have been some attempts 
in the literature to calculate disc models in which these 



processes coexist (e.g. Casse & Ferreira 2000 Ogilvie 



& Livio 2001), these studies have not yet linked the 



generation of radial angular-momentum transport with 
its origin in magnetic effects via MRI-induced turbu- 
lence. In order to devise a method to explore the possi- 
ble joint operation of these two processes, we make use 
of the fact that the MRI is generally suppressed when 
the ratio of the Alfven to sound speed a is of order unity 



(Balbus & Hawley 1991). The instability is expected to 



develop into full turbulence when a 1, provided that 
the field-matter coupling is sufficient to sustain MHD- 
driven processes. On the other hand, wind-driving discs 
require an ordered, relatively strong field (e.g. consis- 
tent with a < 1) for the wind to be launched. Since 
the value of a = _B/[(47rp) 1 / 2 c s ] increases with height 
on account of the strong decline in density away from 
the midplane, this suggests that in principle it is possi- 
ble for both forms of transport to operate at the same 
radius, with the MRI transporting angular momentum 
radially outwards over a finite section of the disc close 
to the midplane, and vertical transport operating at 
higher z. There may also be a magnetically-inactive 
(dead) zone across the disc midplane, where the mag- 
netic diffusivity is large enough to prevent the magnetic 
field from effectively coupling to the fluid and the disc 
shear. 

The above possibility was first explored by |Salmeron| 
(2007) under the assumption of pure ambipolar 



et al. 



diffusivity and constant field-matter coupling. The ob- 
tained, preliminary, results seem to indicate that there 
is little room for the joint operation - at the same ra- 
dius, but at different heights - of these two mechanisms. 
This is because for this to occur the parameter a must 
satisfy two conditions: First, it must be small enough 
near the midplane to allow MRI modes to grow over a 
finite section of the disc. Second, it has to be sufficiently 
large near the surface for the wind to be launched. 
Physically-viable wind solutions also satisfy 2Aa 2 > 1 
(WK93, KSW10). In practice, these constraints seem 
to require comparatively strong magnetic fields at the 
disc midplane (consistent with ao > 0.5), coupled with 
low (and decreasing with z) values of the field-matter 
coupling parameter, for both transport mechanisms to 
operate in the described fashion. These conditions may 
be difficult to satisfy over extended radial locations in 
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real discs. Despite these initial results, I emphasise that 
it is crucial to conduct more detailed modelling, incor- 
porating a realistic diffusivity and field-matter coupling 
profiles, in order to fully explore the relationship be- 
tween radial and vertical angular-momentum transport 
in protoplanetary discs. 



7 Closing remarks 

In this paper I have examined the magnetic activity of 
protostellar accretion discs, the properties of the field- 
matter diffusivity in weakly ionised environments, and 
the two most generally applicable mechanisms for angu- 
lar momentum transport in these objects: MRI-driven 
magnetic turbulence and winds launched centrifugally 
from the disc surfaces. In the first process, the MRI 
redistributes angular momentum radially outwards via 
Maxwell stresses associated with small-scale, disordered 
magnetic fields for which the ratio of the Alfven speed 
to the sound speed is much smaller than unity (a <C I, 
a 'weak' fielcQ. The second mechanism transfers an- 
gular momentum from the disc material to the wind 
through the action of a large scale, ordered field (sat- 
isfying a < 1 already at the disc midplane; a 'strong' 
field). In this latter scenario, magnetic turbulence is 
expected to be suppressed as the MRI-unstable modes 
do not fit within the thickness of the disc (e.g. WK93). 

The low conductivity of the gas, particularly in the 
inner regions of the disc, has to be carefully accounted 
for when modelling these processes. The results indi- 
cate that, despite the weak ionisation, the field is able to 
couple to the gas and shear for fluid conditions thought 
to be satisfied over a wide range of radii in these discs. 
When the Hall diffusivity dominates, however, the mag- 
netic activity depends critically on the magnetic field 
polarity. The likelihood of a joint operation of the two 
angular-momentum transport mechanisms was briefly 
discussed. Further work, incorporating realistic fluid 
conditions, is necessary in order to fully explore this 
important possibility. 

I close with a summary of the main points discussed 
in the paper. 

Magnetic activity 

1. Angular momentum redistribution (hence, accre- 
tion) in protostellar discs is likely to be mediated 
by the action of magnetic fields through the follow- 
ing two processes: MHD-induced turbulence (via the 



8 Note that, in this scenario, the fluctuating component of the 
field may typically be stronger than the mean field (e.g. |Sano| 
|et al. |2004) . 



magnetorotational instability - MRI) or magneto- 
centrifugal jets accelerated from the disc surfaces. 
In protostellar discs, the low ionisation fraction 
makes it imperative to fully account for departures 
to the ideal-MHD limit when evaluating the ability 
of the magnetic field to drive these processes. 
Three diffusivity regimes can be identified in weakly- 
ionised environments: Ambipolar, Hall and Ohm. 
Ambipolar (Ohm) diffusivity dominates in low- 
density (high-density) regions endowed with a strong 
(weak) magnetic field, so that the charges are pri- 
marily tied to the magnetic field (to the neutrals). 
The Hall regime is dominant for intermediate den- 
sities and field strengths, and is characterised by 
a varying degree of coupling of different charged 
species with the magnetic field: Small charges (such 
as electrons) are coupled to the field, whereas more 
massive ones (e.g. ions and dust grains) are tied to 
the neutral component of the fluid. 
In a minimum-mass solar nebula disc, the ambipo- 
lar, Hall and Ohm regimes are dominant in the disc 
interior on radial scales r > 10 AU, ~ 1 — 10 AU, 
and ~ 0.1 — 1 AU, respectively. 
When the Hall diffusivity is important, the mag- 
netic response of the disc is no longer invariant un- 
der a global reversal of the magnetic field polarity 



(e.g. Wardle fc Ng 1999 SKW11) 



Magnetocentrifugal Jets 

6. Angular momentum can be transported vertically, 
via the acceleration of magnetocentrifugal jets from 
the disc surfaces. For this to occur, the magnetic 
field must be relatively strong, so that the ratio of 
the Alfven speed to the sound speed at the midplane 
is of order unity (ao < 1), and the magnetic coupling 
must be sufficiently strong to sustain MHD-driven 
processes. 

7. Winds accelerated centrifugally from the surfaces of 
strongly-magnetised discs exhibit the following lay- 
ered structure, in order of increasing height above 
the disc midplane. 

(a) a quasi-hydrostatic layer, where the bulk of the 
matter is located, matter loses angular momen- 
tum to the field and strong gradients in the fluid 
variables take place; 

(b) a transition layer, where the magnetic energy 
becomes dominant over the fluid thermal energy 
on account of the drop in density away from the 
midplane; and 

(c) a wind region, where the field transfers angu- 
lar momentum back to matter and accelerates it 
centrifugally 
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In the top-most layer, which constitutes the base of 
the outflow, all velocity components [(v — vk4>)/Cs] 
axe positive and increasing with z. In contrast, in 
the first two (e.g. within the disc), the radial veloc- 
ity is negative, the vertical velocity is small and the 
azimuthal velocity (v^) is sub-Keplerian. 

8. We have developed a procedure to explicitly solve for 
the vertical structure of a radially-localised wind- 
driving disc model (WK93, KSW10). The solu- 
tions are integrated numerically between the mid- 
plane and the sonic surface of the flow, and incor- 
porate the strong stratification of the disc (WK93, 
SKW11). Furthermore, by matching the solution 
to the BP82 global self-similar wind models, we are 
able to obtain solutions that continue to accelerate 
to the Alfven surface. 

9. In order to facilitate the matching of the localised 
(wind- driving disc) solution to the global, self-similar 
(large-scale wind) solution, we obtained a library of 
the global solutions of BP82 for a wide range of the 
model parameters. The parameter combinations for 
viable solutions are available in tabular form in the 
electronic version of SKWI1. 

10. We have derived the following parameter con- 
straints, satisfied by physically- viable wind-driving 
disc solutions (WK93, KSW10). 

(a) The flow remains sub-Keplerian below the launch-| 
ing height. 

(b) A wind-launching criterion is satisfied. This im- 
plies the sufficient inclination of the magnetic 
field lines with respect to the rotational axis of 
the disc for centrifugal acceleration to take place. 

(c) The base of the outflow is located above one 
magnetically-reduced density scale- height. 

(d) The rate of heating by Joule dissipation at the 
midplane does not exceed the rate of gravita- 
tional potential energy released at that location. 

11. We find that the above constraints lead to four pa- 
rameter sub-regimes where solutions are found in the 
Hall limit, and to three in the Ohm limit. Our nu- 
merical solutions are in agreement with these con- 
straints. In the regions of parameter space excluded 
by the constraints, wind-driving disc solutions can- 
not be obtained or are unphysical. 

12. For all diffusivity regimes, the ratio of the neutral- 
ion momentum exchange time to the disc orbital 
time is larger than unity (T > 1). This is a fun- 
damental property of the wind solutions considered 
here. 

13. Increasing the relative contribution of the Hall dif- 
fusivity term, while keeping the remaining model 
parameters unchanged, results in a higher location 
of the sonic point, and in a smaller magnetically- 
reduced density scale-height. Correspondingly, the 



density at the sonic point, and the mass outflow rate, 
are both reduced. 

14. When the Hall diffusivity is important, both the 
properties of the solutions and the regions of param- 
eter space where these are found are dependent on 
the magnetic field polarity. In two Hall sub-regimes 
(Cases ii and iv, KSW10; see also Table [T]) wind so- 
lutions exist only for the positive polarity (B z > 0, 
meaning that the field is aligned with the angular 
velocity vector of the disc). In the other two (Cases 
i and Hi), no viable solutions are found for ?7h < 
when T < 2|/3 i0 |- 1 (see Fig. [|]). 

15. When solutions with both field polarities exist, the 
base of the wind and the sonic surface are both lo- 
cated at a higher z in the solutions with B z > (for 
similar values of all the other model parameters). 

Magnetorotational instability 

16. The magnetorotational instability (MPJ) is the 
most promising angular-momentum transport mech- 
anism in discs that are endowed with a relatively 
weak magnetic field, so that the ratio a of the Alfven 
speed to the sound speed is <C 1, provided that the 
coupling between the magnetic field and the gas is 
enough to support MHD-driven processes. 

17. In the weakly- ionised environment typical of the 
inner regions of protostellar discs, the magnetic dif- 
fusivity is strong enough to significantly modify the 
development and properties of the instability. 

18. A local, linear analysis of the instability (WS11) 
shows that for the same value of the Pedersen dif- 
fusivity r]p, the growth rate increases with the Hall 
diffusivity term, from zero when 77h^k/w^ < —2 to 
~ 0.75 Hk, the growth rate associated with ideal- 
MHD conditions, as rfn — > oo. 

19. Hall diffusion modifies the extent of the magnetically-| 
active column by one order of magnitude, increasing 

it (decreasing it) when the magnetic field is parallel 
(antiparallel) to the angular velocity vector of the 
disc (WS11). 

20. MRI-unstable perturbations grow in protostellar 
discs for a wide range of fluid conditions and mag- 
netic field strengths. When dust grains are absent, 
the MRI grows at r = 1, 5 and 10 AU for B < 8 
G, B < 800 mG and B < 250 mG, respectively 
( |Salmeron fc Wardle |[2005| ). 

21. At radial distances within a few AU from the cen- 
tral object, a small (submicron-sized) population of 
dust grains renders the midplane magnetically inac- 
tive. When the magnetic field is aligned with the 
rotational axis of the disc (B z > 0, resulting in 
r/H > 0), this effect is restricted to a finite range of 
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grain sizes. Once the grains have aggregated to ~ 1 
micron or more, the midplane is again unstable to 
the MRI ( |Salmeron fc Wardle |2008| WS11). If the 
field is counter-aligned, however, only a small section 
of the disc is unstable, even if the grains aggregate 
significantly, or are absent (WS11). 
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